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ABSTRACT 

We report on the methodology and first results from the Deep Lens Survey (DLS) transient search. 
We utilize image subtraction on survey data to yield all sources of optical variability down to 24*^ 
magnitude. Images are analyzed immediately after acquisition, at the telescope and in near-real 
time, to allow for foUowup in the case of time-critical events. All classes of transients are posted to 
the web upon detection. Our observing strategy allows sensitivity to variability over several decades 
in timescale. The DLS is the first survey to classify and report all types of photometric and as- 
trometric variability detected, including solar system objects, variable stars, supernovae, and short 
timescale phenomena. Three unusual optical transient events were detected, flaring on thousand- 
second timescales. All three events were seen in the B passband, suggesting blue color indices for the 
phenomena. One event (OT 20020115) is determined to be from a flaring Galactic dwarf star of spec- 
tral type dM4. From the remaining two events, we find an overall rate of 77 = 1.4 events deg~ day" 
on thousand-second timescales, with a 95% confidence limit of 77 < 4.3. One of these events (OT 
20010326) originated from a compact precursor in the field of galaxy cluster Abell 1836, and its 
nature is uncertain. For the second (OT 20030305) we find strong evidence for an extended extra- 
galactic host. A dearth of such events in the R passband yields an upper 95% confidence limit on short 
timescale astronomical variability between 19.5 < Mr < 23.4 of rjR < 5.2 events deg^^ day^^. We 
report also on our ensemble of astrometrically variable objects, as well as an example of photometric 
variability with an undetected precursor. 

Subject headings: gamma rays: bursts — minor planets, asteroids — stars: variables: other — super- 
novae: general — surveys 



1. OPTICAL ASTRONOMICAL VARIABILITY 

Characterization of the variable optical sky is one of 
the new observational frontiers in astrophysics, with vast 
regions of parameter space remaining unexplored. At 
the faint flux levels reached by this optical transient 
search, previous surveys were only able to probe down 
to timescales of hours. An increase in observational sen- 
sitivity at short timescale and low peak flux holds the 
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promise of detection and characterization of rare, violent 
events, as well as new astrophysics. 

The detection of transient optical emission provides 
a window into a range of known astrophysical events, 
from stellar variability and explosions to the mergers of 
compact stellar remnants. Known types of catastrophic 
stellar explosions, such as supernovae and gamma-ray 
bursts (GRBs), produce prompt optical transients decay- 
ing with timescales of hours to months. Some classes of 
GRBs result from the explosion of massive stars - these 
hypernovae are known to produce bright optical flashes 
decaying with hour-long timescales due to emission from 
a reverse shock plowing into ejecta from the explosion. In 
addition, GRBs produce optical afterglows, decaying on 
day to week timescales, resulting from jet-like relativis- 
tic shocks expanding into a circumstellar medium. Even 
more interesting are explosive events yet to be discovered, 
such as mergers among neutron stars and black holes. 
These may have little or no high-energy emission, and 
hence may be discove rable only at longer wavelengths 
ijLi fc Paczvnskil[l998() . Finally, there is the opportunity 
to flnd rare examples of variability, as well as the poten- 
tial for discovering new, unanticipated phenomena. 

One of the primary science goals of this transient 
search is the rate and distribution of short timescale as- 
tronomical variability. However, we catalogue and re- 
port all classes of photometric or astrometric variability 
with equal consideration. A weakness in many surveys 
with targeted science is that serendipitous information 
is discarded as background - a counterexample is the 
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wealth of stellar variability information being gleaned 
from microlensing survey data. The DLS survey geom- 
etry and observation cadence are not optimal for max- 
imizing the ov erall number of detected transients (e.g. 
lNemir off '2003') , being driven instead by weak lensing sci- 
ence requirements. However, given the current lack of 
constraints at short timescales and to significant depth, 
the potential remains for the discovery of new types of 
astronomical variability. 

The search for variability without population bias is 
one of the primary goals of the Large-aperture Synoptic 
Survey Telescope (LSST, ^vson 2002). Efforts such as 
the DLS transient search are a useful and necessary test- 
ing ground for the software design and implementation 
needed to reduce, in real-time, LSST's expected data 
flow of 20 TB nightly. 

2. VARIABILITY SURVEYS 

The field of optical variability surveys has been revo- 
lutionized by the adaption of CCD devices to astronomy, 
and by Moore's progression of computing power. The 
former has yielded wide field imaging systems on large 
aperture telescopes, and the latter ensures their data 
can be reduced in quick order, if not real-time. Such 
advances have lead to modern variability surveys char- 
acterized by successively more extreme combinations of 
depth, cadence, and sky coverage. 

The advent of microlensing surveys advanced by orders 
of magnitude the known number of variable objects, driv- 
ing the forefront of wide-field (0.1-1 deg^) imaging on 
meter-cla ss telescopes (e.g. EROS: Afonso et al. 200i 
MAC HO: lAlcock et ahfc oOOb: MEGAT lAh^eset ani200a 
MOA: iSumi et al.'<2003: OGLE LII.HL iWozniak et all 
l20Mt POINT-AGAPE: Paulin-Hcnriks son et alJ 12003 
SuDerMACHO: »Rest et al.<l2004i|) . The cadences are typ- 
ically of order 1 image per day per field, although vari- 
ous strategies have been employed for sensitivity to short 
events or fine structure in longer events. Such surveys are 
traditionally directed at crowded systems of stars serving 
to backlight a foreground microlensing population. This 
stellar background, however, also serves as a foreground 
that obfuscates sources at cosmological distances. 

Surveys for supernovae at cosmological distances avoid 
such crowded foregrounds, employing similar resources 
and attaining similar sky coverage, but reaching deeper 
than the mic rolensing searches (e.g. High-z: Tonrv et al. 
l2003t SCP: iKnop etal] 1200^ . The original surveys' 
cadence of field reacquisition was typically of order a 
month, short enough to avoid confusion by active galactic 
nuclei but long enough to ensure a handful of new super- 
nova per observing run. These surveys were designed 
only for event discovery, with subsequent lightcurve cov- 
erage provided by narrower field foUowup resources. 
Given the successes of the original supernovae surveys, 
subsequent efforts have been allocated resources for an 
advanced observing cadence allowing both detection and 
self-foUowup, su ch as the IfA Deep ijBarris et alJl200:^ 
and ESSENCE ijSmith et a,lJl2002bD supernova surveys. 
The Nearby Supernova Factory is expecting ~ 100 type 
la supernova per year from their moderate aperture 
1.2-meter) telescope and large (9 deg^) field of view 
Aldering et al . 2002). However, dedicated systems with 
a small field of view (0.01 deg^), similar depth, and a ca- 
dence of several days can lead to a wealth of supernova 



detections (e.g. LOSS: iFilippenko et al.ll200H) . 

Asteroid and Near Earth Object (NEO) searches 
cover of order thousands of degrees per night, primar- 
ily in the ecliptic, and near opposition (e . g. LONEOS: 
'Howell et al.' 1996; LINEAR: iStokes e t al.' 2000; NEAT: 
iPravdo et al. 1999) . Their photometric depths are shal- 
lower than supernova surveys and the cadence more 
rapid, given the requirement of recovering objects in mo- 
tion. The images are typically of poorer seeing and image 
sampling, given less stringent photometric requirements. 

Truly wide-field systems (10 to 100 deg^) incorpo- 
rating small apertures and limiting survey depths of 
12-17*'' magnitudes are able to cover the entirety of 
the sky in a given night, or a given large region 
of sky multiple instances per nigh t (e.g. GROSCE : 
lAkerlof et al. 1993: ROTSE-III : Akerlof et^fl 1200^ 
RAP TOR: IVestrand et all l2002 ASAS-3: |Poimahski| 
12001(1 . Such systems are optimized for the detection and 
followup of fast optical transients, such as optical coun- 
terparts to GRBs. These systems compromise Nyquist 
sampling of the point spread function (PSF) and pho- 
tometric depth for breadth and real-time agility. RAP- 
TOR'S advanced implementation includes a narrow-field 
fovea to optimally study localized t ransients. The North- 
ern Sky Variability Survey fNSVS: 'Wozn iak et alJl200|) 
makes use of ROTSE-I survey data, which sampled the 
entire local sky up to twice nightly and provides variabil- 
ity information on ~ 14 million objects between 8-15.5*'' 
magnitudes. 

Additional variability surveys include the QUEST RR 
Lyrae survey (Vivas et al. 2001), which is to cover 700 
deg2 toV = 2P* magnitude in the Galactic equator, with 
sensitivity to several timescales of variability. The Sloan 
Digital Sky Survey has a sequence of multi-color tempo- 
ral information for over 700 deg^ at multiple timescales 
(Ivezic ct al. 2003). The Palomar-Quest synoptic sky 
survey drift-scans 500 deg^ per night, w ith scans sepa- 
rated by time baselines of days to months l|Mahabal et alJ 
|2003|) . The Deep Ecliptic Survey monitors 13 deg^ per 
night near the ecliptic in a search for Kuiper Belt Objects 
(KBOs), using 4-meter telescopes and 2 exp osures sepa- 
rated by 2-3 hours and a third within a day jMillis et alJ 
| 2002tl. Finally, th e Faint Sky Variability Survey (FSVS: 
iGroot et alJ 120031) is a study of overall optical and as- 
trometric variability at faint magnitudes on a 2.5-meter 
telescope. Point source limiting magnitudes reach to 25 
over a field of ~ 0.1 deg^, and the cadence of field reacqui- 
sitions is staggered for sensitivity at several timescales, 
and thus to multiple phenomena. Of all on-going vari- 
ability surveys, the FSVS is most similar to the DLS 
transient search. 

3. THE DEEP LENS SURVEY 

The Deep Lens Survey (DLS) is a 5-year NOAO survey 
operating on the 4-meter Blanco and Mayall telescopes 
-f- MOSAIC imagers at the Cerro Tololo and Kit t Peak 
observatories, respectively l|Wittman et alJl2002|) . The 
survey is undertaking very deep multicolor imaging of 
five 2° X 2° fields chosen at high Galactic latitude, and 
composed of 9 subfields apiece. Exposure times and fil- 
ters are chosen to reach limiting magnitudes oi B,V, R, z' 
to 29/29/29/28 mag per square arcsecond surface bright- 
ness, with images typically acquired near new moon. At 
these limiting magnitudes, we expect to measure ac- 
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curately the shapes and color-redshifts of ~ 100,000 
galaxies per square degree for weak lensing science. An 
overview of the survey and public data releases are avail- 
able at http : / / dls . bell-labs . com/ The survey is ex- 
pected to conclude observations in March, 2005. 

Full exposure depth requires 20 exposures per pass- 
band in each of our 45 subfields. The images are ob- 
tained, with some exceptions, in sets of 5-pointing dither 
sequences, with offsets of 100" and 200" around the ini- 
tial central pointing to fill gaps between MOSAIC chips 
and enable super-sky-flat construction. Typical expo- 
sure times are 600 seconds in B, V, and z', and 900 sec- 
onds in R, yielding limiting magnitudes per exposure for 
point sources of approximately 24. The weak lensing sci- 
ence goals of the survey require co-adding these images. 
In the co-addition process, however, variable sources are 
clipped out, or averaged away. Here we pursue an or- 
thogonal approach, which instead analyzes differences in 
these data, at the pixel level, to extract this variability. 
This paper summarizes our transient search to date. 

4. DLS TRANSIENT SEARCH 

Starting in December 1999, we began analysis of our 
acquired images in near-real time, which we define as 
within the duration of a dither sequence of observa- 
tions. True real-time reductions ideally occur within the 
timescale it takes to acquire the subsequent observation, 
or, given a target population of events, shorter than the 
timescale of the variability being searched for. Recog- 
nizing variability in near-real time with current CCD 
mosaics requires data analysis of several GB nightly of 
images on-site. To this end, we maintain computers in 
the control rooms of both the Kitt Peak and Cerro- 
Tololo 4-meter telescopes. At Kitt Peak, we have a 
quad-processor 700 MHz machine, and at Cerro-Tololo 
a quad-processor 550 MHz machine. Both systems run 
under the Linux operating environment, and execute the 
custom data reduction pipeline described in Section 

The ensemble of Deep Lens Survey images considered 
for this publication was taken from 14 observing runs be- 
tween 2002 November and 2003 April where the transient 
pipeline was operating efficiently. Figure ^ displays the 
distribution (dotted line) of temporal intervals between 
subfield reacquisitions, integrated over all subfields. The 
solid line shows this same quantity, but only for reac- 
quisitions in the same ffiter as the previous observation. 
Since variability is detected through comparison of ob- 
servations taken in similar filters, the integral under this 
solid subfield-filter histogram up to timescale r defines 
our overall sensitivity to transients at timescale r. The 
integral under the dotted histogram represents our fol- 
lowup capability at a given r. 

Figure ^ indicates our primary sensitivity is to vari- 
ability on thousand-second timescales. In this subset we 
have a total integrated exposure time of 10.2 days dis- 
tributed over 425 S-band, 464 F-band, and 393 i?-band 
images. At 0.6° x 0.6° per image, this yields a total ex- 
posure E of 3.7 deg^-days for sensitivity to lO'^ second 
events down to 24*'* magnitude. We have a factor of 
17 smaller E for inter-night observation intervals at 10^ 
seconds, and also at inter-month observation intervals at 
10^ seconds. A planned re-analysis of all survey data for 
transients will allow for a precise determination of E at 
all available timescales. 



4.1. Pre-pipeline Image Calibration 

Our typical subfield observation comes as a 5- 
image dither pattern controlled using the IRAF^^ script 
mosdither. Standardized information is written in each 
of the images' FITS headers, and is used in the pipeline 
to automatically associate each quintuple of images. 

Basic image calibration is done before each image 
is sent to the transient pipeline. The IRAF task 
mscred . ccdproc is used to correct for cross-talk, apply 
overscan correction and trim the image, and bias and 
flat-field correct the images (we generally take calibra- 
tion sequences in the afternoon to allow for real-time 
reductions). In the case of z'-band data, we need to 
apply a fringing correction, requiring an accumulation of 
images that generally precludes real-time reductions. Fi- 
nally, the images are registered to the World Coordinate 
System (WCS) on the s ky using mscred.m sccmatch with 
the N0A0:USN0-A2 f Monet et al1 IT99Sl catalog. The 
images are then copied to a directory that is continuously 
searched by the transient pipeline image reaper. 

4.2. Automatic Pipeline 

Our transient pipeline uses the OPUS environment 
ijRose et al.lll995ft as a backbone sequencing a series of 
Python and Perl language scripts, each defining a "stage" 
in our pipeline. These stages are defined as follows: 

• IN : copy an image via the SCP protocol from the 
IRAF reduction machine to the transient pipeline 
machine. 

• RI : register and pixel-wise resample an incoming 
image to a fiducial template image of the subfield. 
WCS information is used as a starting point for 
determining the registration coefficients. The as- 
trometric template is the same for all passbands, 
and is constructed from the very first observation 
of a subfield. 

• CI : detect and catalog objects in the registered 
image using SExtractor (iBertin fc Arnouta.1996.) . 
For our detection thresholds, we use SExtractor 
parameters DETECT_THRESH = ANALYSIS.THRESH = 
2.5. 

• HI : difference each dither against t he first point- 
ing, using a modified version of the lAlardI l)200nf l 
algorithm. Modifications include internal methods 
to register and remap pixels based upon WCS in- 
formation, comprehensive discrimination of appro- 
priate regions to determine the image convolution 
kernel, figures of merit to determine the direction 
of convolution, robust methods to ensure complete 
spatial constraint on the variation of the kernel, 
and the propagation of noise through the convolu- 
tion process 

• FI : catalog the difference image and filter object 
detections, inclu ding the calculation of a daptive 
second moments llBernstein fc Jarvi^l2002|) . These 

IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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moments, originally designed for accurate shape 
measurements in our weak lensing pipeline, are 
used here to verify the integrity of each candidate. 
Cuts on the results of this analysis (e.g size criteria 
for cosmic ray rejection) help to reject a significant 
fraction of the noise and false positives, as well as a 
fraction of the actual transients present. This stage 
yields the short-timescale transients in each dither 
sequence, primarily solar system objects, but also 
including stellar and unknown sources of variabil- 
ity. 

To algorithmically cull false detections from our 
object lists, we first construct exclusion zones 
around bright star s, using the N0A0:USN0-A2 
ijMonet et alJll998j) catalog, and along bad CCD 
columns, using bad pixel masks supplied with the 
mscred package. Given a list of candidates SEx- 
tracted from the difference image, we further ana- 
lyze each object in both the original input image 
and the difference image. In the input image, we 
require the object's semi-minor axis, as computed 
from adaptive second moments, to be greater than 
1 pixel (0.27") as a cut against cosmic rays. In 
the difference image, we further require second mo- 
ments greater than 0.9 pixels, uncertainties in the 
orthogonal ellipticity components less than 0.5 pix- 
els, and integrated flux measurements greater than 
counts. We also examine all pixels within 20 pix- 
els of a candidate, and require 50% of all pixels to 
be neither negative nor masked, and 70% of non- 
masked pixels to be positive valued. These addi- 
tional cuts safeguard against bad registrations or 
subtractions, which tend to leave a number of con- 
tiguous pixels above the detection threshold, but 
also regions of systematic dipole residuals. This 
does impact our efficiencies for slowly moving ob- 
jects such as KBOs. By correctly propagating noise 
arrays through the remapping and difference imag- 
ing process, it is possible to compare projected 
noise properties of the image to empirical measure- 
ments, and automatically identify and reject such 
regions. 

At this point, the individual images are co-added to 
make a deeper representation of the subfield and searched 
for longer timescale transients : 

• SI : co-add the dither sequence of images into a 
single deep image of the subfield, providing com- 
plete spatial coverage by filling in MOSAIC chip 
gaps. 

• IT : identify prior templates of the same subfield, 
in the same filter, to use as templates for detecting 
variability on longer timescales. 

• HS : difference the deep image stack constructed in 
stage SI against the prior templates. 

• FS : catalog and filter the deep difference images. 
Longer timescale transients are revealed at this 
stage, primarily supernova and active galactic nu- 
clei (AGN) , and stellar variability. We employ sim- 
ilar object cuts as in stage FI. 



• CL : clean up and archive data files, and inform the 
observers that the sequence of observations is ready 
to be reviewed. 

Our automated cuts typically reduce the number of 
candidates per image from several hundreds to several 
tens. After the last stage of the pipeline, the observer 
must manually categorize the remaining candidates. A 
graphical user interface displays the sequence of differ- 
ence images. Each is presented as a triplet of postage 
stamps centered on the candidate and showing the prior 
template, current image, and difference between them. 
The nature of the transient is determined using a vari- 
ety of factors, which can be divided into temporal and 
contextual information. Temporal discriminants include 
astrometric motion of the transient, indicative of solar 
system origin, as well as the timescale of photometric 
variability for stationary objects, which can help dis- 
criminate between populations of variability with known 
timescales, or to single out an object varying with un- 
expected rapidity. Contextual discriminants include the 
existence of a precursor or host object, and the proximity 
and relative location of the transient event. 

During visual classification of a transient, a set of data 
files are automatically generated locally on disk, after 
cross-checking for spatial coincidence with previous tran- 
sients, or proximity to approximate ephemerides for mov- 
ing objects. An ephemeris is generated for moving ob- 
jects, and for stationary objects, a finding chart. These 
files are copied to our remote web server where a peri- 
odic daemon executes a command to regenerate our web 
page based upon the data files present. Thus, every 15 
minutes our web page is brought into sync with the tran- 
sients detected at the telescope. At the conclusion of each 
run, all moving objects are reported to the Minor Planet 
Center^^ and the most convincing supernova candidates 
to the Central Bureau for Astronomical Telegrams^"*. 

The magnitudes resulting from photometry on differ- 
ence images are denoted by the symbol A4. This rep- 
resents the magnitude of a static object whose bright- 
ness is equal to the measured change in flux A/, M = 
—2.5 log (A/) -|- mo, where toq is the photometric zero- 
point of the template image. This quantity can be re- 
lated to the change in magnitude AM of an object with 
quiescent magnitude M through the relation 

/ (M-A1) \ 

AM = -2.5 log nO^^ + Ij . (1) 

In the limit where the variability comes from a precursor 
fainter than the detection limit, such as for supernovae, 
then A4 represents the magnitude of the transient itself. 
In the following, we report all magnitudes derived from 
difference imaging in units of A4 . 

We are further pursuing algorithms that might auto- 
matically and optimally discriminate between transient 
classes, which will be a necessity for larger scale transient 
searches. One such example is the GENetic Imagery Ex- 
ploration software package fGENIE: lPerkins et al.12 000'1. 
GENIE uses genetic algorithms to assemble a sequence 
of image operators that optimally extract features from 
multi-spectral data, and might be naturally be applied 
instead to our multi-temporal sequences. 
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All archived trans ients are made av ailable on the inter- 
net at http://dls.bell-labs.coin/transients.html 

generally within a few hours of observations. Transients 
are broadly categorized as "moving" (slow or fast objects 
are highlighted), "variable star", "supernova", "agn or 
supernova" (if this cannot be distinguished based on the 
host morphology and transient location), or "unknown", 
a broad category indicating objects whose astrophysical 
nature is uncertain. This includes very rapidly varying 
objects, or transients with no obvious precursor or host. 

5. ASTROMETRICALLY VARIABLE OBJECTS 

Roughly 75% of the transients detected in the DLS 
are astrometrically variable solar system objects. Four 
of the five DLS fields are within 15° of the ecliptic, and 
in these fields we generally find of order ten moving ob- 
jects per subfield visit. The total number of moving ob- 
jects detected through spring 2003 is 3651, with 900 de- 
tected in B, 1508 detected in V, and 1487 detected in R 
(some were detected in multiple filters). Our detection 
efficiency for moving objects has not been modeled. 

Because of our observing pattern, most of the moving 
objects are observed on a single night and never recov- 
ered. A main belt asteroid moving at ~30" hr~^ near 
opposition leaves a subfield within a few days. As a re- 
sult, only about one-third of our moving objects have 
been observed on multiple nights, and 333 of these have 
received Minor Planet Center provisional designations, 
which are listed on our website. Similarly, we are able 
to measure colors for only a fairly small fraction of these 
objects. Therefore, we do not attempt a thorough anal- 
ysis of orbital parameters, albedos, and sizes. Rather, 
we give an overview of the dataset and then detail a few 
especially interesting objects. 

Figure 121 shows the magnitude distribution of all mov- 
ing objects discovered, separately for B, V, and R. Note 
that the magnitude range is displaced from the magni- 
tude range in which our efficiency for stationary point 
sources is nonzero (Figure ^J. This is because the vast 
majority of moving objects are significantly trailed dur- 
ing our 600 or 900 second exposures. An asteroid at the 
bright end can therefore be much brighter than a point 
source before saturation sets in, and at the faint end must 
also be brighter to rise above the detection threshold. 

Figure |3| shows the moving objects' velocity distribu- 
tion in ecliptic coordinates, highlighting a pair of Kuiper 
Belt objects (KBOs) detailed below. Figure 0] provides 
tentative classifications into asteroid families (main belt, 
Trojans, Centaurs, NEOs, etc.) ba sed upon ecl i ptic v e- 
locities, adopting the boundaries of llveziis "eTaLl |200l|) . 

5.1. Kuiper Belt Objects 

The 2003 April CTIO run yielded two Kuiper Belt Ob- 
jects (KBOs). Our pipeline does not include an orbit cal- 
culator, so KBOs are not automatically flagged. These 
two examples were selected manually on the basis of their 
small (< 3" hr~^) angular motions near opposition. 

Transient 207 from the 2003 April CTIO run was dis- 
covered on 2003 April 1 (UT), moving at 2.9" hr~^ near 
opposition, wit h a median R magnitude of 21.5. Using 
the method of IBer nstein & Khushalani (2000), we find 
a seminiajor axis, a, of 41.0 AU, eccentricity e of 0.022, 
and inclination i of 20.1°, making this a classical KBO. 
Assuming an albedo of 0.1, the diameter of this object 



is ~475 km, making it one of the larger KBOs known. 
Transient 587 from the same run was discovered on 2003 
March 31 (UT) at a _B magnitude of 23.3, moving at 2.6" 
hr~^, also near opposition. Our orbit calculation yields 
a = 42.8, e = 0.178, and i = 1.1°, making it a likely 
scattered KBO. For an albedo of 0.1, this object is ~410 
km in diameter. Our observations provide orbital arcs of 
5 and 4 days for these KBOs, respectively. 

5.2. Near-Earth Objects 

The detection efficiency for fast movers declines with 
velocity, as their reflected light is trailed over larger area. 
The DLS, with its long exposures, is therefore far from 
optimal for detecting Near-Earth Objects (NEOs). Nev- 
ertheless, we have detected a sizable sample, as shown in 
Figure ^ Here we highlight one that was moving with 
an unusual velocity. Transient 1363 from the 2003 March 
KPNO run. This object is the fastest solar system object 
seen during the transient search, moving at 80" hr~^. 
The box near the top of Figure 0] indicates this object's 
remote position in the ecliptic velocity diagram. 

6. PHOTOMETRICALLY VARIABLE OBJECTS 

We focus here on those transients that have been 
noted to flare on thousand second timescales and are 
detected in at least two images, as well as on longer- 
term variability with no detectable originating host or 
precursor. However, our transient catalog also includes 
a selection of several hundred variable stars, as well as 
over 100 supernova candidates, 18 of which pass the 
requirements for recognition as defined by the Central 
Bureau for Astronomical Tel egrams. These include : 
2000bj ( Kirkman etldl l2003 ). 2000fq fWittma n et"al 
12000), 2002aj, 2002ak, 2002al, 2002am (Be cker et al. 
|20.Q2c), 2002ax, 2 002ay, 2002az 2002b a, 2002bb, 2002bc, 
2002bd, 2002be iBecker et al.l l20q2aD . 2003bx, 2003by, 
2003bz, 2003ca (jBecker et al.ll2003j) . 

6.1. Image Calibration 

We have re-run the difference imaging algorithm for 
this paper, using 2k x 2k (8.6' x 8.6') arrays centered 
on the selected optical transients (OTs). We have gen- 
erally been able to construct deeper, higher signal-to- 
noise (S/N) template images than those available to us 
at the time of detection. For photometric calibration, 
we need only calibrate the zero points of these tem- 
plate images, as the difference imaging convolution nor- 
malizes each input image to the template's photometric 
scale. The B, V, and R templates are normalized to the 
iLandolti ()1992|) system, with additional errors (typically 
5%) added in quadrature to account for instrumental dif- 
ferences between the KPNO, CTIO, and Landoh (199^ 
systems. Al l z' tem plates are initially calibrated to the 
iSmith et alJ l)2002a|) AB system, and presented here in 
the Vega system by adopting an offset of 0.54 magni- 
tudes {z'y^gg^ = z'y^g — 0.54). The transformation be- 
tween photometric systems depends on the instrumental 
filter response and on Vega's spectral energy distribution, 
adding an additional 0.08 magnitude uncertainty to our 
z' = z'y^g^ measurements. 

The template brightnesses of the precursor, or host, 
obje cts were determined usin g MAG_BEST from SExtrac- 
tor l)Bertin fc Arnout^ Il996(l . except in the case of 
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OT 20010326 (Section EXH below). For this precur- 
sor, which appears unresolved, we first used the IRAF 
noao . digiphot . daophot . psf package to determine the 
local PSF from nearby isolated stars, and subtracted 
these stars and the transient precursor from the image 
using noao . digiphot . daophot . allstar. 

For those template images where we are not able to de- 
tect the OT precursor/host, we define the point source 
detection limit in accordance with the NOAO Archive 
definition of photometric depth (T. Lauer, private com- 
munication), as 



mi 



TOO - 2.5 log(1.2 W CTsky n), 



(2) 



where mg is the magnitude of one ADU, W is the seeing 
in pixels, fisky is the dispersion of the image around its 
modal sky value, and we choose n = 3 cr as our detection 
limit. We determine detection limits in the difi^erence 
images by measuring directly our ability to recover input 
steUar PSFs. 

For the difference imaging, we choose the best -fit con- 
volution kernel that varies spatially to order one across 
our 2k X 2k subimage, and a sky background with a 
similar degree of variation. We fit these spatial vari- 
ations over a 20 X 20 grid, using a convolution kernel 
with a half-width of 10 pixels. The function of the con- 
volution kernel is to degrade the higher quality image 
to match the PSF of the lesser quality image. Objects 
are detected in each grid element to constrain locally 
the convolution kernel through a pixel-by-pixel compar- 
ison between images. A global fit is next performed, 
using each element as a constraint on spatial variation 
of the kernel. Several sigma-clipping iterations ensure 
that sources of variability, such as variable stars or as- 
teroids, are not used as constraints for this kernel. If such 
an element is rejected from the global fit, a replacement 
object within this element is chosen to ensure maximal 
spatial constraint. Application of the final convolution 
kernel to the entire image and pixel-by-pixel subtrac- 
tion yields the difference between the two images, with 
a resulting point spread function of the lesser quality 
image, and photometric normalization to the template 
image. The transient fluxes were determined using the 
IRAF noao . digiphot . apphot .phot package. The rel- 
ative transient magnitudes are typically determined to 
better than ~ 4%, with the remaining uncertainty due 
to absolute calibration of the template zero-point. 

6.2. Optical Transients Without Hosts 

We have detected several optical transients with no ob- 
vious hosts, and which arc seen to vary over the timescale 
of months (e.g. iBeckcr 2003a), similar to supernovae 
and GRB afterglows. One of the more lumi nous exam- 
ples of this class of objects is OT 20020112 (jBecker et al.l 
l2002bt) . originally reported as Transient 139 from our 
2002 January run Information on this event is sum- 
marized in Tabled including the DLS subfield in which 
it was detected (F4pl3), and position on the sky in J2000 
coordinates. We include limits on a precursor or host for 
this event, which is not detected to deeper than 27*'' mag- 

-th . 



assum ing a Ry — 3.1 extinction curve ijSchlegel et alJ 
119981) . 

OT 20020112 was detected in observations taken 2002 
January 12. Previous subfield observations were taken 
2001 March 27, so the age of the transient at detec- 
tion is poorly constrained. Information on this tran- 
sient was quickly posted to the Variable Star NET work 
(VSNET^^), prompting a sequence of foUowup observa- 
tions from the 1.3m McGraw-HiU telescope (MDM). The 
lightcurve for this transient is shown in Figure \^ and 
represents a daily averaged lightcurve of all observations 
from the DLS and MDM. The event was re-imaged 35 
days after detection as part of DLS survey observations, 
and had faded and reddened significantly. Detailed pho- 
tometry for this event is listed in Table |2 The mag- 
nitudes listed are in differential flux magnitudes M - 
however, given the non-detection of a precursor or host, 
this also represents the traditional magnitude M of this 
transient. 

OT 20020112 was also observed spectroscopically with 
the Las Campanas Observatory 6.5-m Baade telescope 
(-1- dual imager/spectrograph LDSS2) on 2002 January 
18. The spectrum in Figure El is characterized by a blue 
continuum with no obvious broad features, and marginal 
evidence for emission lines Ha at z = 0.038 and [Nil] 
6583 at z = 0.039. The continuum features and lack 
of short timescale variability are consistent with a Type 
II SN caught close to maximum brightness. However at 
z = 0.038, the observed brightness (F - S - 22.1) is 
3 magnitudes dimmer than a typical Type II SN near 
maximum. The blue continuum makes a reddened SN 
explanation unlikel y. Examinatio n of the FIRST 1.4 GHz 
radio catalog (, White et alJll99'^ yields a limit of 0.94 
mjy/beam at this position, with no object closer than 
2.4^ A search of the ROSAT source catalog fVoees et ajj 
120001) . and a more comprehensive search of archival X- 
ray and gamma-ray data through HEASARC^^ yields no 
archival sources within 10'. 

Thus the nature of this object remains unknown, a sit- 
uation compounded by the lack of an obvious host galaxy 
to i? > 27.6. We note the observed B — V color change 
during the event is inconsistent with a spectral energy 
distribution of the form o<i , such as that expected 
for G RB afterglows (e.g. iSari et al.lll9M iNakar et alJ 
l2002j) . However, evolution of the index (i could lead to 
the observed behavior. Future OT searches will require 
facilities dedicated to foUowup of such events, to asso- 
ciate or distinguish them from the expected supernova 
and GRB afterglow populations. 

6.3. Rapid Variability Optical Transients 

Over the course of the survey we have also detected 
three optical transients on thousand-second timescales 
that were present in multiple images, and whose pre- 
cursors were not immediately recognized as stellar in all 
available passbands. In all cases, we were able to identify 
a precursor or host object after the fact. These are des- 
ignated OT 20010326, OT 20020115, and OT 20030305. 



nitudes in the B, V, and R passbands, and ~ 25"" magm 
tude in z'. Tabled also lists limits on the precursor/host 
brightness after taking into account Galactic reddening. 



|http: //vsnet .kusastro .kyoto-u. ac ■ jp/vsnet/index.html| 

^"^ A service of the Laboratory for High Energy Astrophysics 
(LHEA) at NASA/ GSFC and the High Energy Astrophysics 
Division of the Smithsonian Astrophysical Observatory (SAO), 
|http: //dls .bell-labs ■ com/transients/ Jan- 2002/transl39 .html http : //heasarc . gsf c . nasa . gov/db-perl/W3Browse/w3br owse ■ plj 
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Tabled lists the observed characteristics of the precur- 
sors or hosts for these transients in their quiescent state. 
Tabled also includes magnitudes corrected for Galactic 
reddening for OT 20010326 and OT 20030305, in the case 
that they lie outside the Galactic dust layer. 

When putting our subsequent results in a cosmolog- 
ical context, we assume a WMAP cos mology of Ha ~ 
71, f^Af = 0.27, Qa = 0.73 llSpergel et al.ll2(¥)3i) . 

6.3.1. OT 20010326 : Mar 2001, Transient 52 

OT 20010326 was detected 2001 March 26.2 as the 52"'' 
transient during that particular run Photometry for 
this event lightcurve, including the constraining obser- 
vations immediately preceding and following the event, 
is listed in Table |21 After verification, an alert was dis- 
persed via e-mail on Mar 26.4 requesting immediate fol- 
lowup observations. 

The event lightcurve is shown in Figure [3 and is char- 
acterized by a null detection in B, followed by a i3-band 
peak and immediate fall-off. Given the range of allowed 
turn-on times, from 100s before the start of the detec- 
tion observation (to account for read-out time between 
images) to Is before the shutter closes in the detection 
image, we can limit a power-law index for the flux decay 
/oc^-" of 0.80<a< 1.2. 

Radio observations of this OT were undertaken at the 
VLA on 2001 March 30. Limits on radio flux at 8.5 GHz 
are —0.1 ±0.3 mJy. Subsequent observations of the event 
revealed a host galaxy or precursor object visible in the 
V , R, and z' passbands, and a limit of _B > 26.4. 

This event occurred in the field of galaxy cluster Abell 
1836, at a redshift of z = 0.037. Figure [7| also shows an 
8.6' X 8.6' i?-band template image centered on the tran- 
sient event. The precursor/host is displayed at higher 
resolution in the upper-right corner as the dim object 
3.5" north-west of the brighter star. A single archival 
HST image of cluster galaxy PKS 1358-11 also covers the 
location of the observed transient. This 500s integration 
in the F606W filter was obt ained Jan 25, 1995 a s part 
of a nearby AGN survey bv IMalkan et all (|1998). The 
precursor appears unresolved in this image. We compare 
this precursor position with the location of two back- 
ground galaxies, and over a baseline of 7.0 years find 
proper motion limits of 0.004 ± 0.004"yr~^. 

It is particularly difficult to photometer and calibrate 
this OT precursor/host, as it occurs 3.5" from a bright 
star that is saturated, or nearly so, in all of our images. 
There is also a very strong background gradient from the 
nearby elliptical galaxy PKS 1358-11 and spiral galaxy 
LCRS B135905.8-112006. If the OT host is a member 
of Abell 1836, it lies only 49 kpc in projection from the 
core of PKS 1358-11, and 53 kpc in projection from the 
center of LCRS B135905.8-112006. The HST resolution 
limit at the distance of Abell 1836 is around 70 parsecs, 
which does not immediately preclude a globular cluster 
host. However, at a distance modulus of 36.0, the host 
absolute magnitude of My = —11.7 would be consider- 
ably brighte r than the globular cluster systems in our 
own Galaxy ijVan D en Bcrgh 2003). addition, the ob- 
served V—R = 1.1 color of the object makes it unlikely to 
be an unresolved globular cluster system (K.A.G. Olsen, 
private communication) . 

^^ |http://dls .bell-labs ■com/traiisients/Mar-2001/trans52.html| 



6.3.2. OT 20020115 : Jan 2002, Transient 337 

Our second fast optical transient, OT 20020115, was 
initially reported as Transient 337 from the run ongo- 
ing 2002 January 15.3^^. The lightcurve, whose data are 
presented in Table |21 and displayed in Figure |H1 is very 
similar to that of OT 20010326, including B-band de- 
tection, but it reaches a brighter A4b, and decays much 
more rapidly. The flux decay power-law index a is con- 
strained to 1.5 < a < 2.4. 

This event was quickly reported as a potential GRB op- 
tical count erpart to the GRB Circular Network (GCN, 
IBecker etl a.1. 2002d). We obtained spectra of the source 
soon thereafter - Figure |51 shows 3 x 10 minutes of ex- 
posure on the source with the Magellan 1 telescope -I- 
LDSS2 dual image/grism spectrograph, taken 2002 Jan- 
uary 18.3. The appearance of zero-redshift emission 
and absorption features in these spectra are evidence 
for Galactic or igin, and the GCN alert was amended 
ijClocchiatti et al. 2002 1. 

Subsequent observations in R and z' revealed a very 
bright, red precursor object, consistently saturated in R- 
band images. We also recover the precursor object in the 
B and V bands. With a quiescent magnitude of i? = 21.3 
and measured peak A^s — 20.7, Equationnindicates the 
precursor varied by more than 1 magnitude, averaged 
over the duration of our 600s exposure. 

The analysis of the photometric observations and ob- 
ject spectrum (Figure 0) shows OT 20020115 has the 
colors and spectrum of a dM4 post-UV Ceti type flare 
star. The spectrum shows the clear presence of late- 
type M-dwarf spectral indicators: molecular bands of 
TiO, CaOH (5530-60A and centered at 6250A), plus 
CaH at 6385A, in addition to weak Ha and H/3 in emis- 
sion. With a V magnitude of 19.67 ± 0.05, and colors 
B-V = 1.58 ±0.09 andV-z' = 3.39 ±0.10, these data 
are consistent with an M dwarf of spectral type dM4, i.e. 
a late type flare star in the quiescent state. Assuming a 
typical My = 12, OT 20020115 would he at a distance 
of approximately 350 pc, and 250 pc above the Galactic 
plane, well within the scale height for disk population M 
dwarfs. 

6.3.3. OT 20030305 : Mar 2003, Transient 153 

The most dynamic transient, OT 20030305, was de- 
tected as Transient 153 from our 2003 March run^°. This 
is the only event detected in the 1^-band, and we imme- 
diately acquired a sequence of S-band observations after 
the y-band series. The event lightcurve is shown in Fig- 
ure ^1 Also displayed are the series of images of this 
event. The top row represents the precursor/host of the 
event in its quiescent state, in the R, I, and z'-bands, 
from left to right. Included in each panel is the magni- 
tude of this precursor/host in its quiescent state, except 
for the / band where we have no calibration data. This 
precursor/host is not visible in B or V, and appears unre- 
solved in the z' images but measurably elliptical in both 
/ and R. The subsequent panels, reading from left to 
right and top to bottom, show the emergence of the OT 
in our time series of images, as well as OT magnitudes 
and time elapsed from the detection observation. All im- 
ages are separated by approximately 700s except for the 



http: //dls .bell-labs ■ com/transients/ Jan- 2002/trans337 .ht ml| 
http: //dls .bell- labs ■ com/trajisient s/Mar-2003/traiis 153 .html | 
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final image, which was obtained more than 2 days after 
the event. The event w as also announced via the lAU 
Circulars l|Beckeill2003'b|) . 

The OT lightcurve shows more complex behavior than 
the previous events, and cannot be described as simple 
power-law flux decay. The second temporal peak indi- 
cates either a second round of brightening or spectral 
evolution of the flare. 

As noted above, the precursor/host object appears re- 
solved and elliptical in the R and / passbands. The 
object's position angles in R and /, derived with the 
adapt ive second moment method of iBernstein fc .Tarvial 
l)2002D . agree to within 0.09°, indicating an extended ob- 
ject. We have characterized the expected PSF at the 
location of the host object by first identifying candi- 
date stars, based on their position in the size-magnitude 
plane, in the i?-band template image. We then fit for 
the spatial variation of the moments of these candidates, 
clipping at 3a to reject interloping small galaxies. A 
test comparing the shape of the host with the predicted 
shape of the PSF at its position yields — 11-4 for 3 
degrees of freedom, for a confldence level of 99% that a 
point source would not have the measured shape of the 
precursor object. We consider the above as strong evi- 
dence for a resolved, extragalactic host for OT 20030305. 

6.4. Discussion 

A primary discriminant between Galactic stellar 
events, such as OT 20020115, and extragalactic events 
is the measured shape of the precursor or host object. 
OT 20010326 is shown to be unresolved in HST imaging, 
implying Galactic origin, although potentially originat- 
ing from a compact extragalactic host. OT 20030305 has 
an apparent resolved host that is inconsistent with the 
i?-band PSF interpolated at its position. This argues 
strongly for an extragalactic origin. 

It is suggestive that all of our OTs are detected in the B 
passband. We consider various stellar and extragalactic 
possibilities for the origin of these events in turn. 

6.4.1. Galactic Flare Stars 

Since we are certain OT 20020115 is from a flaring 
dwarf star, we investigate the possibility that one, or 
even both, of the other OTs are similar in nature. We 
examine the peak of the disk dwarf luminosity function, 
which is occupied by dM4 - dM6 {My = 12 - 16) stars 
(jReid et al.lll995j) . This region is occupied by the most 
active UV Ceti-type flare stars, with event timescales of 
order min utes and flare B — V colors of a pproximately 
0.0 to 0.3 l)Kunkellll975j) . iGurzadianI l|T980() further char- 
acterizes the "Type I" outburst lightcurves expected of 
UV Ceti-type stars, with rise times of seconds to a few 
minutes, and decay times of minutes to about one hour, 
qualitatively simil ar to our observed even ts. 

From Figure 2 of lBoeshaar et alJ lj2003j) . and FigurelTTI 
here, the precursor R — z' color of OT 20010326 im- 
plies a spectral type of dM4 were it a flare star. With 
B — R > 3.1, this object lies la from the color-defined 
locus of disk main sequence stars. However, using a V- 
band dM4 distance modulus of between 11.5 and 12.5 
magnitudes, this object would be approximately 2.0-3.2 
kpc distant, and 1.5-2.4 kpc above the Galactic plane. 
This instead suggests association with the sdM halo sub- 
dwarf population, consistent with our proper motion lim- 



its. However, these stars are not known to flare. The sit- 
uation is similar for OT 20030305, which is far too faint 
&tV> 27.1 mag to belong to the disk population. As- 
suming an i? — z' derived spectral type of approximately 
dM6, its distance modulus yields d > 1.7 kpc. 

For a typical DLS field at b(II) = 45°, we expect to 
find approximately 100 - 150 dM4 - dM6 disk stars per 
subfield out to their Galactic scale height of 350 pes, 
which is verified through direct examination of several 
subfields. Since fiare events do not follow Poisson statis- 
tics, and show large scatter in total energy, peak light 
and duration, it is difficult to calculate the exact number 
expected in each sequence of B exposures. In addition, 
selection efi^ects in most previous studies have biased dis- 
covery towards stars of greatest flare visibility. 

Until recently, one of the only ways to estimate the 
expected rate of flare st a r acti vity was to extrapolate 
Equation 5 from iKunkell l)1973(l . This relation was de- 
rived using data for a selected few of the most active 
flare stars in the immediate solar neighborhood. Scal- 
ing Kunkel's relation to our B passband, and assuming 
that 5% of all dM4-6 stars exhibit strong flaring activity 
(i.e. AS of 1 mag or greater), we may place an up- 
per limit of approximately 0.08 flare events expected per 
subfield in each 600 sec exposure. Given the extreme se- 
lection bias in Kunkel's sample, it is not unexpected to 
discover that this estimate {N < 34) is much different 
than what we find (1 < < 3). We also examine the 
stellar flare survey data of lFresneau et al. (200lj). Utiliz- 
ing astrographic plates covering 520 deg^ at low Galac- 
tic latitude to apparent B magnitude of 10-14, they find 
8% of their stars show flare events of > 0.4 mag over 
20-30 minutes. If half of these stars are M dwarfs, of 
which approximately 1% are types dM4-6 based on the 
volume of the sample limited by the luminosity of those 
spectral types, we expect an upper limit of 0.01 flare of 
> 0.4 mag per 100 stars in each 600 sec i3-band expo- 
sure. Having analyzed 425 accumulated i?-band images 
for this publication, we expect an upper limit of 6.4 flar- 
ing events > 0.4 magnitudes. The distribution function 
of the amplitudes of flaring events in these stars is not 
well known, but in general the brighter flares make up a 
small er fraction of the total number of events ijGurzadianI 
I1980D . Thus the discovery of only one certain flare event 
of AB > 1 mag is consistent with this upper limit. Given 
the lack of unbiased flare star samples down to our limit- 
ing magnitudes, a direct comparison with other surveys 
requires th e analysis of more modern variability datasets 
(e.g FSVS: [Groot et al.ll2f)f)l . 

6.4.2. Active Galactic Nuclei (AGN) 

IVanden Berk et all lj2002|) have reported a highly lumi- 
nous OT event from the Sloan Digital S ky Survey, with 
an un derlying host galaxy at z = 0.385. iGal-Yam et alJ 
(j2002j) have determined, using in part archival optical 
and radio data, that the event is likely due to a radio- 
loud A GN. Heeding the suggestions of IGal-Yam et alJ 
ll2?)0l . we search archival data sources for observations 
of the fields of OT 20010326 and OT 20030305. We in- 
clude in these searches the FIRST 1.4 GHz radio cat- 
alog ( Wh ite et all |1997[ ). 1.4 GHz NVSS radio cata- 
log (Condon et al. 199^, and 4.8 GHz PMN catalogs 
ifGriflith fc Wright.19931 

The only survey that covers the field of OT 20010326 
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is the PMN tropical survey, which detects no sources 
closer than 1.7'. In combination with our VLA detec- 
tion limits 4 days after the event, this argues against a 
radio variable source. Null detections are also reported 
for OT 20030305 in the FIRST, NVSS, and PMN equa- 
torial catalogs, with an explicit limit from FIRST of 0.94 
mJy/beam. Overall, these null detections seem to reject 
radio-loud AGN as sources for this variability. 

We have also searched the archival X~ray and gamma- 
ray catalogs through HEASARC, and while OT 20010326 
returns several matches corresponding to components of 
Abell 1836, no matches are found within 2'. A search 
around OT 20030305 finds no matches within 10'. This 
weighs against a variety of accretion scenarios, including 
AGN and other Galactic X-ray sources. 

6.4.3. Gamma Ray Bursts (GRBs) 

There is an expected population of optical transients 
associated with the observed rate of gamma ray bursts. 
However, the connection between the two event rates is 
very uncertain, given the variety of means to draw optical 
emission from hydrodynamic evolution of a GRB event. 
These means include "orphan" GRBs, where a highly col- 
limated GRB points at large enough angle away from the 
observer that gamma ray emission is avoided but optical 
afterglow radiation is not (Rhoads 1997); and "dirty" 
or "failed" fireballs whose ejecta comprise a significant 
amount of bar yons and/or have a small Lorcntz factor 
(< 100, e.g. iDermer et all IT9'99. Huans ct al. 2002). 
Distinguishing between these possibilities requires fine 
constraints on the time decay slope a, radio observations, 
and/or multiwavelength obs ervations at th e time of the 
so-called GRB "jet break" l|R,hoadsll200l . Given the 
our poor temporal resolution, and lack of simultaneous 
multi-wavelength coverage, we are not uniquely sensitive 
to the micro-physics expected to drive the early evo- 
lution of GRB lightcurves, and whose resolution might 
distinguish between the above possibilities. Within the 
resolution and timescale of our observations, all GRB 
variants effectively yield the same characteristic isolated 
fading object. 

The thousand second timescale of the observed events 
is of similar duration t o the early lig htcurve break ob- 
served in GRB 021211 l|Li et alJl200l . After this break, 
the power law flux decay index a was seen to decrease 
from 1.8 to 0.8 - before this time, the transient dimmed 
by more than 2.5 magnitudes. If our already faint, 
rapidly declining optical transients represent a similar 
early phase of gamma-ray dark GRB lightcurves, their 
subsequent evolution would be below our detection limit. 

We also note that the spectral energy distribution 
(SED) of the GRB afterglow compone nt is modeled as 
a synchrotron emission spectrum (e.g. iSari et 311119981 
iGranot et al.1 ^0001. The expected spectral shape scales 
like Fi, (X . Thus it is surprising that our distri- 
bution of i?:l/:i?-detected transients is 2:1:0 (where OT 
20030305 contributes to both the B and V event rate) 
if these events are drawn from a population of events 
characterized by a synchrotron SED. 

Recent submillimeter and radio observations of local- 
ized GRB host galaxies show they are bluer than se- 
lected galaxies at similar redshifts, indicati ve of star for- 
mation or relatively low du st content (e.g. iBerger et al.l 
120031 iLe Floc'h et alJl200^ . Host galaxies for localized 



X-ray Flashes (XRFs) also appear to exhibit the same 
characteristics (Bl oom et al. 200i|). However, our puta- 
tive hosts are considerably redder than these samples of 
GRB or XRF-selected host galaxies, and are unlikely 
to be drawn from the same population. The R — z' 
colors of our OT hosts are, however, loosely consistent 
with the observed population of Extremely Red Objects 
(EROs), which comprise old or dusty star-forming galax- 
ies. If EROs represent highly obscured starburst activity 
at moderate redshift (z ~ 1), they contribute a signifi- 
cant fraction of th e overall star formation at that epoch 
ijSmail et alJl200^ . Assuming our population of OTs 
traces star formation, we expect some association with 
EROs. The absence of extremely red GRB host galaxies 
suggests detection of GRBs is biased against dusty galax- 
ies. If our OTs arise from dust enshrouded galaxies, this 
bias would seem lesser for optically detected events. 

7. EVENT RATE 

We map our detection of these objects into a quanti- 
tative statement describing the rate of short timescale 
astrophysical variability. We derive general constraints 
without bias towards a particular OT model. 

7.1. Efficiency 

To account for inefficiencies in our pipeline, we model 
our efficiency £ at recovering transients throughout a 
range in brightness. We calibrate the efficiency via Monte 
Carlo runs in which we add point source transients. We 
define an efhciency run as the generation of 32 random 
amplitude point sources, with the integrated flux chosen 
as a uniform deviate in the exponent oi f — 10^-^-^-^, 
which are placed within a dither sequence of survey data. 
Positions are randomly selected within the limits of the 
primary exposure, and are placed at the same astronom- 
ical position within each of the subsequent 4 dithers. We 
do not bias our results by assuming a temporal shape for 
the variability. Instead, we generate an overall efficiency 
for recovering a given amount of input flux, averaged over 
the spatial coverage of our dither sequence. 

We include in our analysis B, V, and i?-band survey 
data (z'-band data are infrequently reduced in real-time 
due to fringing complications). We account for varia- 
tions in observing conditions by choosing 19 night -filter- 
subfield combinations of survey data to submit for effi- 
ciency analysis. In total, we initiated 8028 efficiency runs 
yielding 256896 total efficiency points, approximately 
85000 per passband. The magnitud es of these transients 
are calibrated to the lLandoTO l)1992() system. Image zero- 
point offsets were calculated from directly calibrated im- 
ages of each subfield. One subfield was calibrated using 
publicly av ailable NO AO Deep Wide-Field Survey data 
ijJannuzi fc Dev . 1999) . 

After placing additional flux in each image, we pass 
them through our difference imaging pipeline. The num- 
ber of input objects recovered is tallied, searching both 
the unfiltered and filtered object catalogs generated at 
the FI stage in our pipeline. This yields an unfiltered 
detection efficiency {dotted histogram in Figure IT^ as 
well as the efficiency of passing our cuts (solid histogram 
in Figure lT^ . The difference between the two histograms 
demonstrates the rejection of true positives in our filter- 
ing process, yielding an overall decrease in efficiency. We 
use the solid histogram in the subsequent event rate anal- 
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ysis. We note that approximately 3% of our sky coverage 
lands in gaps between MOSAIC chips, and the dithering 
procedure shifts 10% and 20% of the template images' 
common sky area off of the subsequent dithers. Thus 
our maximum possible efficiency is ~ 82%. Our overall 
OT detection efficiency of £ ~ 63% indicates our filtering 
efficiency is ^ 77%. We expect substantial improvements 
in the context of a future automated classification envi- 
ronment. 

Figure 1121 shows our efficiency is fairly constant be- 
tween the saturation limit on the bright end, and our 
detection limits on the faint end. We thus expect little 
bias against faint transient detection in this given range. 
On the bright end, there appears to be non-zero £ for 
objects brighter than saturation. These are systematic 
artifacts that result from the detection of unsaturated 
and unmasked wings of saturated efficiency objects, in 
close enough proximity to the input object itself to war- 
rant a positional match. In the following analysis, we 
set £^ = brighter than the average saturation limits of 
Mb = 18.6, A^y = 18.8 and Mr = 19.5. 

7.2. General Constraints on Variability 

To average over unknown inefficiencies in the human 
element of our transient pipeline, we require that a tran- 
sient be confirmed in a subsequent image, which we con- 
sider a human observer 100% efficient at classifying as 
real. Thus our limits on short timescale variability cor- 
respond to twice our typical exposure time, plus 100s of 
readout between images. Our efficiencies in Figure 1121 
are relatively constant between the bright and faint lim- 
its, which we can approximate by a constant efficiency 
{£) without loss of specificity. However, due to the dif- 
fering {£) per passband, differences in typical exposure 
times, and in number of events detected, we also quote 
event rates for each passband individually. The observed 
rate for each passband, averaged over all combinations 
of detectable variability amplitudes, is 

77 ~ — ^ — events deg~^ day~^ (3) 

{£f E & y y I 

where N is the number of observed events, E the ap- 
propriate exposure, and {£) is the averaged efficiency, 
squared due to the requirement of two detections. We 
note that OT 20030305 contributes to both the V and B 
event rate, as it was discovered in both sets of images. 

For the B passband with 3 detected transients, our 
overall rate for variability on 1300s timescales, between 
18.6 < Mb < 23.8, is t? = 6.5 events deg-^ day"\ 
We recognize at least one of these events as Galactic in 
origin. Having detected no more than two cosmological 
events, Poisson statistics exclude at the 95% level any 
OT model that predicts a mean number of detectable 
cosmological events Nb > 6.3, implying 773 < 14 events 
deg-2 day~\ Table El lists, for a given number of con- 
sidered events, experimental rates and limits on rj. We 
have detected zero short timescale events in the R pass- 
band, constraining overall 1900s astronomical variabil- 
ity between 19.5 < Mr < 23.4 to 7/95% < 5.2 events 
deg~^ day^^. 

In our total summed exposure of 3.7 deg^-days from all 
passbands, assuming {£) = 0.63 and 4 events (3 unique 
events), we find an overall rate of short timescale astro- 
nomical variability of r] — 2.7 (2.0) events deg~^ day~^. 



with 7795% < 6.3 (5.3). These are the first general con- 
straints on short timescale variability at such depths. 

8. CONCLUSIONS 

We have reported on the structure of and first results 
from our wide-field image subtraction pipeline. An im- 
portant characterization of our transient survey is the ex- 
posure i? at a given timescale and to a given depth. The 
DLS transient search is primarily sensitive to ^ 1000s 
variability from 19*'* to deeper than 23'"'* magnitudes in 
B, V, and R. Within this envelope of sensitivity, we have 
detected three short timescale optical transient events. 

OT 20010326 occurred in the field of galaxy cluster 
Abell 1836. One archival HST image of the cluster in- 
cludes this region, and the transient precursor is present 
and unresolved. This indicates a compact precursor, stel- 
lar in nature if it resides in our Galaxy. The colors of this 
object are consistent at the la level with those expected 
of Galactic dwarf stars fFigure [TT|l . whose flaring activ- 
ity presents a known background. However, the object 
would be too far out of the Galactic plane to belong to 
the disk population of dwarf stars, and halo subdwarfs 
are not known to flare. It is also possible the precursor 
resides in, or behind, the Abell cluster, but overall its na- 
ture remains uncertain. OT 20020115 is identifled as a 
Galactic M dwarf of spectral type dM4, exhibiting classi- 
cal flare star activity. Finally, the host for OT 20030305 
appears consistently elliptical in the R and / passbands, 
and inconsistent with the i?-band stellar PSF, which we 
consider a strong argument for a resolved host, extra- 
galactic in nature. We find OT 20030305 the strongest 
candidate so far for optically detected, short timescale 
cosmological variability. 

The precursor or host objects for OTs 20010326 and 
20030305 are definitively redder than galaxies that are 
known to host GRBs and XRFs. If our events are cos- 
mological in nature, this suggests that optical and high 
energy events arise from different mechanisms, or, given 
the dearth of GRBs from dusty star-forming galaxies, 
a stronger bias against GRBs and XRFs from a dusty 
environment. If our lightcurves evolve analogous to the 
prompt stage of GRBs emission, the rapid decline cou- 
pled with their intrinsic faintness would make them diffi- 
cult to monitor beyond several hours. Finally, we em- 
phasize the diversity of variable objects in the stellar 
menagerie, and thus it is not straightforward to rule out 
Galactic stars as the sources for our OTs. Spectroscopic 
information will ultimately help to clarify their nature, 
and such foUowup observations are planned. 

Our search has also yielded SN-like events that appear 
to have no host galaxy to significant (> 27) limiting mag- 
nitude - OT 20020112 is one such example. In addition, 
our catalog of SN candidates, classified as supernovae 
primarily due to their proximity to a host galaxy, might 
also contain sources with unusual temporal evolution. 

Overall, the DLS transient search is well suited to ex- 
plore the parameter space of OTs with small energy bud- 
gets, a population that could plausibly have escaped de- 
tection by gamma-ray and X-ray satellite missions. This 
raises the possibility that the phenomena detected repre- 
sent a new class of astronomical variability. Coordinated 
photometric foUowup of future optical transients is abso- 
lutely necessary to reveal fully the diversity of variability 
at faint optical magnitudes. A primary goal of future 



11 



variability surveys must be to enable the photometric 
and spectroscopic foUowup of detected events through 
timely release of information and ease of access to avail- 
able data. 

The wealth of information that can be gleaned from 
real-time synoptic transient science, only a subset of 
which is covered in this paper, strengthens the science 
case for the expansion of deep, wide astronomical sur- 
veys into the short timescale regime. Ex pected future 
surveys like the LSST and PAN-STARRS ijKaiser et all 
|2p02) will survey thousands of square degrees per night, 
and will benefit greatly from modern development in this 
field. 
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Fig. 1. — Distribution of sampling intervals between subfield observations using DLS transient survey data. The dotted histogram 
represents this quantity regardless of the filter of observation, whereas the solid histogram represents the same interval restricted to 
identical subfiold-filter comlainations. 
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Fig. 2. — Magnitude distribution of moving objects detected in the DLS. Detections in the B, V, and R passbands are represented by 
dotted, dashed, and sohd lines, respectively. Some objects were detected in more than one filter, but most were not. 




Fig. 3. — The moving object velocity distributions in ecliptic coordinates. Two Kuiper Belt objects detailed in Section Is.ll are marked 
with a red box. They clearly stand out from the main belt distribution, and neighboring points are also likely KBOs. Overall, KBOs 
do not stand out in the distribution. 
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Fig. 5. — Lightcurve OT 20020112 in differential flux magnitudes M in the B (square datapoints), V (triangle), and R (circle) passbands. 
The data were obtained at the CTIO 4-m Blanco telescope and the KPNO 1.3m MDM telescope. All data from a given MJD are averaged 
together after being placed on similar magnitude systems during the difference imaging process. Lightcurve information is presented in 
Table 121 There is no apparent host for this variability in images obtained prior to and subsequent to the observed event. 
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Fig. 6. — Spectrum of OT 20020112 obtained 6 days after detection. The phase of the event is uncertain. Weak emission features indicate 
a possible redshift of z = 0.038. However, we are unable to detect a host to R > 27.6. This transient represents a class of objects that vary 
on timescales of tens of days, much like supernovae or potentially optical GRB orphans, but have no apparent host galaxy. 




Fig. 7. — B and i?-band lightcurve (square and circular symbols, respectively) in differential flux magnitudes A4 (see Equation^, and 
8.6' X 8.6' ij-band template image, of OT 20010326. The latter shows its position in the field of Abell 1836, at z = 0.036. The precursor is 
to the immediate north-west of the (saturated) brighter star, where north is up and east is to the left. Radio source and elliptical galaxy 
PKS 1358-11 is to the south, and spiral galaxy LCRS B135905.8-112006 is to the north-east. 



18 



20 



OT 20020115 



"T 1 1 T 



~i 1 1 r 



"T 1 1 T 



22 



24 



26 



_l I I L 



J I I L 



_l I L 



1000 



1000 
t - MJD 52289.30443631 (sec) 



2000 



Fig. 8. — _B-band temporal lightcurve of OT 20020115 in differential flux magnitudes M. 
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Fig. 9. — Post— event spectrum of OT 20020115. The spectrum was obtained with the Magellan 1 telescope 3 days after the transient 
event depicted in Figure|8] The clear emission in Ha and absorption by TiO bands indicates Galactic origin. Overall, the spectral features 
and composite colors are consistent with a late type (dM4) flare star. 
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Fig. 10. — Lightcurve of OT 20030305, in differential flux magnitudes M, in the B and V passbands (square and triangular datapoints, 
respectively). The image sequence corresponding to these data points is shown in 40" x 40" panels to the right. The top 3 images show the 
host, or precursor, of the transient in R, I, and z' (we have no calibration data for the /-band image). The images below show the time 
evolution of the event, from the constraining image 700s before, to the first image 2 days after where it was undetected in V. The time of 
each observation relative to the detection of the event is indicated, as well as the measured brightness of the transient, where these numbers 
were derived from photometry on the difference images. Our first detection of this transient landed on the boundary of a MOSAIC chip, 
the dark region in the above detection image. In all images, north is up and east is to the left. 



20 




12 3 4 

R-z' 



Fig. 1 1. — B — R vs. R — z' diagram for the precursor /host objects of OT 20010326 and OT 20030305. The open squares represent the 
IPickled llT99^ stellar library of energy distributions, with solid squares indicating the approximate positions of dM3-6 dwarf stars. The 
tracks indicate K-corrected colors for elliptical (circle), irregular (triangle), and Sb (diamond) galaxies with increasing rcdshift (steps of 
z = 1 are indicated by filled symbols, starting with the large open symbol at z = 0) based on galaxy spectral templates from Coleman ct aT] 
119801) . These are shown to indicate regions where contamination of stellar precursor identifications by compact high rcdshift galaxies is 
expected. 
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Fig. 12. — Transient point— source detection efficiency as a function of magnitude in the difference image (M), for the B, V, and R 
passbands. The dotted histograms are our raw point source detection efficiencies, while the solid histograms represent our efficiencies 
after cuts. The bright end cutoff represents the level at which the objects become saturated, with additional efficiency coming from the 
recovery of the wings of the saturated object. In our efficiency analysis, we explicity set our efficiency for objects brighter than saturation 
{Mb < 18.6, Mv ^ 18-8 a^nd Mn < 19.5) equal to zero. The dim cutoff represents the 2.5 a detection limit of single images in our 
transient search. 
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TABLE 1 

Unusual Transient Event Summary 



Event 






DEC ( 12000"! 




y 


R 




z' 


OT 20020112 


F4pl3 


10:49:32.8 


-04:31:45 


> 27.3 


> 27.2 


> 27.7 




> 24.9 










> 27.2 


> 27.0 


> 27.6 




> 24.8 


OT 20010326 


F5p31 


14:01:42.2 


-11:35:15 


> 26.4 


24.48 ±0.10 


23.30 ± 


12 


21.16 ±0.11 










> 26.2 


24.27 ±0.10 


23.13 ± 


12 


21.06 ±0.11 


OT 20020115 


F4p23 


10:48:56.1 


-05:00:41 


21.25 ±0.07 


19.67 ±0.05 


< 19.5 




16.28 ±0.09 


OT 20030305 


F4p31 


10:53:45.8 


-05:37:44 


> 27.2 


> 27.1 


24.58 ± 


10 


21.38 ±0.13 










> 27.0 


> 27.0 


24.50 ± 


10 


21.34 ±0.13 



Note. — Summary of selected optical transients (OTs) from the DLS transient survey. OT 20020112 varied on the timescalc 
of months, and we fail to detect a host or precursor object to the reported hmiting magnitudes. The other three events varied 
on thousand— second timescales, and were resolved in multiple survey images. Magnitudes are Vega-based and represent any 
detected precursor or host in its quiescent state. A second set of magnitudes takes into account^ihlegel et al. I 199^ corrections 
for Galactic reddening, assuming a J?v — 3.1 extinction curve, except for OT 20020115 which is known to be Galactic in origin. 
Point source limits from Equation [2] are quoted for null detections. 



TABLE 2 

Optical Transient Event Photometry 



Event 


MJD 


Passband 


M 


Mc 


OT 20020112 


52285.3 


V 


22.24 ±0.11 


22.11 ±0.11 




52286.2 


B 


22.24 ±0.06 


22.07 ±0.06 




52287.2 


R 


22.10 ±0.08 


22.00 ±0.08 




52287.4 


V 


22.17 ±0.15 


22.04 ±0.15 




52288.2 


R 


22.23 ±0.04 


22.13 ±0.04 




52288.2 


B 


22.32 ±0.08 


22.15 ±0.08 




52288.3 


V 


22.18 ±0.11 


22.05 ±0.11 




52290.2 


B 


22.42 ±0.07 


22.25 ±0.07 




52290.3 


V 


22.31 ±0.11 


22.18 ±0.11 




52319.2 


V 


24.07 ±0.11 


23.94 ±0.11 




52321.3 


B 


25.46 ±0.14 


25.29 ±0.14 


OT 20010326 


51994.2315949 


B 


> 24.8 


> 24.4 




51994.2397245 


B 


22.76 ±0.06 


22.50 ±0.06 




51994.2478576 


B 


23.70 ±0.07 


23.44 ±0.07 




51994.2646898 


R 


> 24.7 


> 24.5 


OT 20020115 


52289.2963074 


B 


> 25.0 






52289.3044363 


B 


20.73 ±0.07 






52289.3125651 


B 


23.83 ±0.07 






52290.2594999 


B 


> 24.8 




OT 20030305 


52703.1732501 


V 


> 23.8 


> 23.7 




52703.1813791 


V 


20.99 ±0.12 


20.89 ±0.12 




52703.1895094 


V 


22.22 ±0.09 


22.12 ±0.09 




52703.1995017 


B 


22.50 ±0.08 


22.37 ±0.08 




52703.2076305 


B 


21.57 ±0.06 


21.44 ±0.06 




52703.2157589 


B 


21.49 ±0.06 


21.36 ±0.06 




52703.2249089 


B 


21.70 ±0.06 


21.57 ±0.06 




52703.2330383 


B 


21.77 ±0.07 


21.64 ±0.07 




52705.3304513 


V 


> 23.9 


> 23.8 



Note. — Differential flux magnitudes for each optical transient, in the 
Vega system. Also listed arc the Modified Julian Day (MJD) and filter for 
each observation. OT 20020112 varied on a month-long timescalc, and the 
photometry listed represents the MJD— averaged brightness. The constraining 
observations preceding and following the event are also listed for the short 
timescalc transients. Limits on A4 are determined by adding input PSFs to 
the difference images. Error bars include systematic calibration errors added 
in quadra ture with the experim ental uncertainties in the measurement. A4c 
represents lSchlegel et alj il998f) corrected magnitudes, assuming a source out- 
side the Galactic extinction layer, except for OT 20020115 which is known to 
be Galactic in origin. 
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TABLE 3 

Event rate 



Filter 


Timescale 


Exposure E (deg^-days) 






{£) 


Nobs 


Nmax,95% 




Vmax,95% 


r> 


ioUUS 


1.1 


10. 


Zo.o 


u.oo 


Q 
O 


7 ft 


0.0 


1 7 
1 / 














2 


6.3 


4.3 


14 














1 


4.7 


2.2 


10 

















3.0 


0.0 


6.5 


V 


1300s 


1.2 


18.8 


23.3 


0.63 


1 


4.7 


2.1 


9.9 

















3.0 


0.0 


6.3 


R 


1900s 


1.5 


19.5 


23.4 


0.62 





3.0 


0.0 


5.2 


B+V+R 




3.7 






0.63 


4 


9.2 


2.7 


6.3 














3 


7.8 


2.0 


5.3 














2 


6.3 


1.4 


4.3 














1 


4.7 


0.7 


3.2 

















3.0 


0.0 


2.0 



Note. — For each filter, we present our sensitivity to transients for timcscalcs at twice the typical exposure times plus 
100s of camera readout. We list the magnitude ranges over which our efficiency of recovering a transient is approximated 
Ijy a constant {£) . Assuming we have detected < A'' < A^j.^^ events per filter, we also list the maximum number of 
events allowed by Poisson statistics at the 95% confidence level, for any model of temporal variability. We derive the 
corresponding event rate r} in events deg~^ day~^. 



